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I shortly review the present status of the theoretical estimates of e' / e . 
I consider a few aspects of the theoretical calculations which may be rel- 
evant in understanding the present experimental results. I discuss the 
role of higher order chiral corrections and in general of non-factorizable 
contributions for the explanation of the AI =1/2 selection rule in kaon 
decays and e'/e . Lacking reliable lattice calculations, the 1/N expansion 
and phenomenological approaches may help in understanding correlations 
among theoretical effects and the experimental data. The same dynamics 
which underlies the CP conserving selection rule appears to drive e'/e in 
the range of the recent experimental measurements. 
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1 Introduction 



The results announced during the last year by the KTeV |1] and NA48 |2| collab- 
orations have marked a great experimental achievement, establishing 35 years after 
the discovery of CP violation in the neutral kaon system || the existence of a much 
smaller violation acting directly in the decays. 

While the Standard Model (SM) of strong and electroweak interactions provides 
an economical and elegant understanding of indirect (e) and direct (V) CP violation 
in term of a single phase, the detailed calculation of the size of these effects implies 
mastering strong interactions at a scale where perturbative methods break down. In 
addition, CP violation in K — > titt decays is the result of a destructive interference 
between two sets of contributions, which may inflate up to an order of magnitude the 
uncertainties on the hadronic matrix elements of the effective four-quark operators. 
All that makes predicting e'/e a complex and subtle theoretical challenge [ffl. 

The status of the theoretical predictions and experimental data available before 
the KTeV announcement in February 1999 is summarized in Fig. [IJ. 

The gray horizontal bands show the one-sigma average of the old (early 90's) 
NA31 (CERN) and E731 (Fermilab) results. The vertical lines show the ranges of the 
published theoretical predictions (before February 1999), identified with the cities 
where most of the group members reside. The range of the naive Vacuum Saturation 
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Figure 1: The 1-a results of the NA31 and E731 Collaborations (early 90's) are shown by 
the gray horizontal bands. The old Miinchen, Roma and Trieste theoretical predictions for 
e'/e are depicted by the vertical bars with their central values. For comparison, the VSA 
estimate is shown using two renormalization schemes. 
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Figure 2: The latest theoretical calculations of e'/e are compared with the combined 1-c 
average of the NA31, E731, KTeV and NA48 results (e'/e = 19.2 ± 4.6 x 1(T 4 ), depicted 
by the gray horizontal band (the error is inflated according to the Particle Data Group 
procedure when averaging over data with substantially different central values). 



Approximation (VSA) is shown for comparison. 

The experimental and theoretical scenarios have changed substantially after the 
first KTeV data and the subsequent NA48 results. Fig. |2| shows the present experimen- 
tal world average for e'/e compared with the revised or new theoretical calculations 
that appeared during the last year. 

Notwithstanding the complexity of the problem, all theoretical calculations show 
a remarkable overall agreement, most of them pointing to a non- vanishing positive 
effect in the SM (which is by itself far from trivial). 

On the other hand, if we focus our attention on the central values, many of the 
predictions prefer the 10~ 4 regime, whereas only a few of them stand above 10~ 3 . Is 
this just a "noise" in the theoretical calculations? 

The answer is no. Without entering the details of the various estimates, it is 
possible to explain most of the abovementioned difference in terms of a single effect: 
the different size of the hadronic matrix element of the gluonic penguin Qq obtained 
in the various approaches. 

While some of the calculations, as the early Miinchen and Rome predictions, 
assume for (ttti\Qq\K) values in the neighboroud of the leading 1/N result (naive 
factorization), other approaches, first of which the Trieste and Dortmund calculations 
and more recently the Lund and Valencia analyses, find a substantial enhancement of 
this matrix element with respect to the simple factorization result. The bulk of such 
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an effect is actually a global enhancement of the 1 = components of the K — > txtx 
amplitudes, which affects both current- current and penguin operators, and it can be 
at least partially understood in terms of chiral dynamics (final-state interactions). 

2 Final State Interactions 

As a matter of fact, one should in general expect an enhancement of e'/e with re- 
spect to the naive VSA due to final-state interactions (FSI). As Fermi first argued [|], 
in potential scattering the isospin 1 = two-body states feel an attractive interaction, 
of a sign opposite to that of the 1 = 2 components thus affecting the size of the cor- 
responding amplitudes. This feature is at the root of the enhancement of the 1 = 
amplitude over the 1 = 2 one and of the corresponding enhancement of e'/e beyond 
factorization. 

The question is how to make of a qualitative statement a quantitative one. A 
dispersive analysis of the K — > -nit amplitudes has been recently presented in ref. [|]. 
The Omnes-Mushkelishvili representation |7j 

Mis + ie) = Pis) exp [ - f° ^ ds') (1) 

\1T J Ami S' — S — IE ) 

is used in order to resum FSI effects from the knowledge of the irn rescattering phase 
5(s) in the elastic regime (s < 1 GeV 2 ). P(s) is a polinomial function of s which is 
related to the factorized amplitude. A solution of the above dispersive relation for 
the A 0t 2 amplitudes can be written as 

Ms) = A' (s - ml) TZj{s) e"'W, (2) 

where A' is the derivative of the amplitude at the subtraction point s = m\. The 
coefficient 1Z represents the rescaling effect related to the FSI. By replacing A'j with 
the value given by LO chiral perturbation theory, Pich and Pallante found 7^.(m|) .2 — 
1.4, 0.9 thus confirming via the resummation of the leading chiral logs related to FSI 
the enhancement of the 1 = amplitudes, together with a mild depletion of the 1 = 2 
components. 

The numerical significance of these results has been questioned M on the basis 
that the precise size of the effect depends on boundary conditions of the factorized 
amplitude which are not unambiguously known, due to higher order chiral corrections. 
As a matter of fact, while the choice of a low subtraction scale minimizes the effect 
of momentum dependent chiral corrections the result of ref. |J cannot account for 
polinomial corrections due to contact interaction terms whose size is unknown (and 
renormalization-scheme dependent) . 
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The analysis of ref. || shows non-perturbatively the presence of a potentially large 
departure from factorization which affects the I = K — > tttt matrix elements. Nev- 
ertheless, the question whether the FSI rescaling of the factorized isospin amplitudes 
leads by itself to a satisfactory calculation of e'/e remains open. 



3 CP conserving versus CP violating amplitudes 

Given the possibility that common systematic uncertainties may affect the cal- 
culation of e'/e and the A J = 1/2 rule (see for instance the present difficulties in 
calculating on the lattice the "penguin contractions" for CP violating as well as for 
CP conserving amplitudes p)) a convincing calculation of e'/e must involve at the 
same time a reliable explanation of the AI = 1/2 selection rule, which is still missing. 
FSI effects alone are not enough to account for the large ratio of the 7 = 0,2 CP con- 
serving amplitudes. Other sources of large non-factorizable corrections are needed, 
which may affect the determination of e'/e as well. 

The AI = 1/2 selection rule in K — > tttt decays is known since 45 years |]TD| and 



it states the experimental evidence that kaons are 400 times more likely to decay in 
the 1 = two-pion state than in the 1 = 2 component. This rule is not justified by 
any general symmetry consideration and, although it is common understanding that 
its explanation must be rooted in the dynamics of strong interactions, there is no up 
to date derivation of this effect from first principle QCD. 

The ratio of I = 2 over 1 = amplitudes appears directly in the definition of 
e'/e : 
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(3) 



As a consequence, a self-consistent calculation of e' / e must also address the determi- 
nation of the CP conserving amplitudes. 

The way we approach the calculation of the hadronic K — > tttt transitions in gauge 
theories is provided by the Operator Product Expansion which allows us to write the 
relevant amplitudes in terms of the hadronic matrix elements of effective AS = 1 
four quark operators (at a scale /x) and of the corresponding Wilson coefficients, 
which encode the information about the dynamical degrees of freedom heavier than 
the chosen renormalization scale: 

G 

h as =i = Ks H [^0) + r y*0*)] Q*0) • ( 4 ) 

The entries of the 3x3 Cabibbo-Kobayashi-Maskawa matrix describe the flavour 
mixing in the SM and r = —VtdV t */V u dV* s . For ji < m c (q = u,d,s), the relevant 
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quark operators are: 



Ql = Oa«,s)v-A ( U P d a)y-A 

Q 2 = (su) Y _ A (ud) y _ A 
Qs,5 = Nv-aE,(?9)v t a 

Qa,6 = (s a dp) v _ A J2q(Qi3 ( la)wT^ 
Ql$ = § (sd) y _ A J^q &q (w)y±A 

Qs.io = | is a dp) y _ A J2 g e q {q f3 q a )y±A 



Current- Current 



Gluon "penguins" 



Electroweak "penguins" 



(5) 



Current-current operators are induced by tree-level W-exchange whereas the so-called 
penguin (and "box") diagrams are generated via an electroweak loop. Only the lat- 
ter "feel" all three quark families via the virtual quark exchange and are therefore 
sensitive to the weak CP phase. Current-current operators control instead the CP con- 
serving transitions. This fact suggests already that the connection between e'/e and 
the AI = 1/2 rule is by no means a straightforward one. 

Using the effective AS = 1 quark Hamiltonian we can write e'/e as 
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(7) 



and (Qi) = (inr\Qi\K) . The rescattering phases <5o,2 can be extracted from elastic 
7i — 7i scattering data JTT]] and are such that cos<5o ~ 0.8 and cos52 — 1- Given that 
the phase of e, 9 e , is approximately 7r/4, as well as 5o — 62, <t> = | + 62 — So — £ turns 
out to be consistent with zero. 

Two key ingredients appear in eq. || 



1. The isospin breaking 7r° — 77 — 77' mixing, parametrized by which is esti- 
mated to give a positive correction to the A 2 amplitude of about 15-35%. The 
complete inclusion of NLO chiral corrections to the tt — 77 — 77' mixing fi2| , |l"3"| and 
of additional isospin breaking effects (AI = 5/2 Ul~4, 151,161) in the extraction of 
the isospin amplitudes may sizeably affect the determination of e'/ e . Although 
a partial cancellation of the new terms in e'/e reduces their numerical impact, 
we must await for further analyses in order to confidently assess their relevance. 

2. The combination of CKM elements ImA t = lm(V t *V t d), which affects directly 
the size of e'/ e and the range of the uncertainty. The determination of Im X t 
depends on B-physics constraints and on e |L7|. In turn, the fit of e depends on 
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the theoretical determination of Bk, the K — K hadronic parameter, which 
should be consistently determined within every analysis. The theoretical un- 
certainty on Bk affects subtantially the final uncertainty on ImA(. A better 
determination of the unitarity triangle is expected from the B-factories and 
the hadronic colliders. In K-physics, the decay — > -n^vv gives the cleanest 
"theoretical" determination of ImA t , albeit representing a great experimental 
challenge. 



4 Summary of theory results 

A satisfactory approach to the calculation of e' / e should comply with the following 
requirements: 

A: A consistent definition of renormalized operators leading to the correct scheme 
and scale matching with the short-distance perturbative analysis. 

B: A self-contained calculation of all relevant hadronic matrix elements (including 
B K ). 

C: A simultaneous explanation of the AJ =1/2 selection rule and e'/e . 

None of the existing calculations satisfies all previous requirements. I summarize 
very briefly the various attempts to calculate e' / e which have appeared so far leading 
to the estimates shown in Figs. [I] and |2|. 

• VSA: A simple and naive approach to the problem is the VSA, which is based 
on two drastic assumptions: the factorization of the four quark operators in 
products of currents or densities and the saturation of the completeness of the 
intermediate states by the vacuum. As an example: 

{ti + 7i-\Q 6 \K } = 2 (n-\u l5 d\0}(7T + \su\K ) -2 (n + n-\dd\0}(0\s-f 5 d\K°) 

+ 2 [(0|ss|0) - (0|deZ|0)] (ti + ti- \s l5 d\K°) (8) 

The VSA does not allow for a consistent matching of the scale and scheme 
dependence of the Wilson coefficients (the HV and NDR results are shown in 
Fig. |l|) and it carries potentially large systematic uncertainties f|]. It is best 
used for LO estimates. 

• Taipei's: Generalized factorization represents an attempt to parametrize the 
hadronic matrix elements in the framework of factorization without a-priori 



assumptions [18]. Phenomenological parameters are introduced to account for 



non-factorizable effects. Experimental data are used in order to extract as much 
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information as possible on the non-factorizable parameters. This approach has 
been applied to the K — > tttc amplitudes in ref. The effective Wilson co- 



efficients, which include the perturbative QCD running of the quark operators, 
are matched to the factorized matrix elements at the scale fip which is arbi- 
trarily chosen in the perturbative regime. A residual scale dependence remains 
in the penguin matrix elements via the quark masses. The analysis shows that 
in order to reproduce the AI = 1/2 rule and e'/e sizable non-factorizable 
contributions are required both in the current-current and the penguin matrix 
elements. However, some assumptions on the phenomenological parameters 
and ad hoc subtractions of scheme-dependent terms in the Wilson coefficients 
make the numerical results questionable. The quoted error does not include any 
short- distance uncertainty. 

Munchen's: In the Miinchen approach (phenomenological 1/N) some of the 
matrix elements are obtained by fitting the AI =1/2 rule at fi = m c = 1.3 
GeV. On the other hand, the relevant gluonic and electroweak penguin (Qq) 
and (Qs)2 remain undetermined and are taken around their leading 1/N values 
(which implies a scheme dependent result). In Fig. [2] the HV (left) and NDR 
(right) results are shown pt)| . The dark range represents the result of gaussian 



treatment of the input parameters compared to flat scannning (complete range). 
Dortmund's: In the recent years the Dortmund group has revived and improved 



the approach of Bardeen, Buras and Gerard |Tl|] based on the 1 /N expansion. 
Chiral loops are regularized via a cutoff and the amplitudes are arranged in 
a p 2n /N expansion. A particular attention has been given to the matching 
procedure between the scale dependence of the chiral loops and that arising 
from the short-distance analysis |p2]] . The renormalization-scheme dependence 



remains and it is included in the final uncertainty. The AI = 1/2 rule is 
reproduced, but the presence of the quadratic cutoff induces a matching scale 
instability (which is very large for Bk)- The NLO corrections to (Qq) induce a 
substantial enhancement of the matrix element (right bar in Fig. |2|) compared 
to the leading order result (left bar). The darker ranges correspond to central 
values of m s , Q v +r)', ImAi and Aqcd- 

Dubna's: In the Nambu, Jona-Lasinio (NJL) modelling of QCD |23[ the Dubna 
group El has calculated e'/e including chiral loops up to 0{p & ) and the effects 



of scalar, vector and axial- vector resonances. Chiral loops are regularized via the 
heat-kernel method, which leaves unsolved the problem of the renormalization- 
scheme dependence. A phenomenological fit of the AI = 1/2 rule implies 
deviations up to a factor two on the calculated (Qe). The reduced (dark) range 
in Fig. corresponds to taking the central values of the NLO chiral couplings 
and varying the short- distance parameters. 
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Trieste's: In the approach of the Trieste group, based on the Chiral Quark 
Model (xQM) p5] , all hadronic matrix elements are computed up to 0(p 4 ) in 
the chiral expansion in terms of the three model parameters: the constituent 
quark mass, the quark condensate and the gluon condensate. These parameters 



are phenomenologically fixed by fitting the AI = 1/2 rule [26]. This step is 
crucial in order to make the model predictive, since there is no a-priori argument 
for the consistency of the matching procedure (dimensional regularization and 
minimal subtraction are used in the effective chiral theory). As a matter of 
fact, all computed observables turn out to be very weakly dependent on the 
scale (and the renormalization scheme) in a few hundred MeV range around 
the matching scale, which is taken to be 0.8 GeV as a compromise between 
the ranges of validity of perturbation theory and chiral expansion. The 1 = 
matrix elements are strongly enhanced by non-factorizable chiral corrections 
and drive e'/e in the 10 -3 regime. The dark (light) ranges in Fig. [| correspond 
to Gaussian (flat) scan of the input parameters. The bar on the left represents 



the result of ref. (g8j which updates the 1997 calculation p7| . That on the right 



is a new estimate El, similarly based on the xQM hadronic matrix elements, 
in which however e'/e is computed by including the explicit computation of 
e in the ratio as opposed to the usual procedure of taking its value from the 
experiments. This approach has the advantage of being independent from the 
determination of the CKM parameters Im A( and of showing more directly the 
dependence on the long-distance parameter Bx as determined within the model. 

Roma's: Lattice regularization of QCD is the consistent approach to the prob- 
lem. On the other hand, there are presently important numerical and theoretical 
limitations, like the quenching approximation and the implementation of chiral 
symmetry, which may substantially affect the calculation of the weak matrix 
elements. In addition, chiral perturbation theory is needed in order to obtain 
K — > 7T7T amplitudes from the computed K — > tt transitions. As summarized 
in ref. lattice cannot provide us at present with reliable calculations of 



the 1 = penguin operators relevant to e'/e , as well as of the 1 = compo- 
nents of the hadronic matrix elements of the current-current operators (penguin 
contractions), which are relevant to the AI = 1/2 rule. This is due to large 
renormalization uncertainties, partly related to the breaking of chiral symmetry 
on the lattice. In the recent Roma re-evaluation of e'/e (Qq) is taken at the 



VSA value with a 100% uncertainty |3(| . The result is therefore scheme depen- 



dent (the HV and NDR results are shown in Fig. ^j. The dark (light) ranges 
correspond to Gaussian (flat) scan of the input parameters. 



Montpellier's: The analysis in ref. [HTJ is based on QCD Sum Rules and uses 



recent data on the r hadronic total decay rates. The value of the Qs matrix 
element thus found is substantially larger than the leading 1/N result. At the 
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same time, the matrix element of the Q§ gluonic penguin, computed assuming 
scalar meson dominance, is found in agreement with leading order 1/N. The 
combined effect is a strong cancellation between electroweak and gluonic pen- 
guins which leads to a vanishingly small e'/e . Various sources of uncertainties 
in the calculation and the comparison with other analyses are discussed in [EHI] . 



Lund's: The AI = 1/2 rule and Bk have been studied in the NJL framework 



and 1/N expansion by Bijnens and Prades |32| showing an impressive scale 
stability when including vector and axial-vector resonances. The same authors 
have recently produced a calculation of e'/e at the NLO in 1/N |33|]. The 
calculation is done in the chiral limit and it is eventually corrected by esti- 
mating the largest SU(3) breaking effects. Particular attention is devoted to 
the matching between long- and short- distance components by use of the X- 
boson method p^ , p5| . The couplings of the X-bosons are computed within the 



ENJL model which improves the high-energy behavior. The A J = 1/2 rule is 
reproduced and the computed amplitudes show a satisfactory renormalization 
-scale and -scheme stability. A sizeable enhancement of the Q§ matrix element 
is found which brings the central value of e'/e at the level of 3 x 10~ 3 . 

• Valencia's: The standard model estimate given by Pallante and Pich is obtained 
by applying the FSI correction factors obtained using a dispersive analysis a 
la Omnes-Mushkelishvili JFj to the leading (factorized) 1/N amplitudes. The 
detailed numerical outcome has been questioned on the basis of ambiguities 
related to the choice of the subtraction point at which the factorized amplitude 
is taken ||. Large corrections may also be induced by unknown local terms 
which are unaccounted for by the dispersive resummation of the leading chiral 
logs. Nevertheless, the analysis of ref. || confirms the crucial role of higher order 
chiral corrections for e'/e , even though FSI effects alone leave the problem of 
reproducing the A J = 1/2 selection rule open. 

Other attempts to reproduce the measured e'/e using the linear cr-model, which 
include the effect of a scalar resonance with m a ~ 900 MeV, obtain the needed 



enhancement of (Qq) [j36[ . However, it is not possible to reproduce simultaneously 
the experimental values of e'/ e and of the CP conserving K — > tctc amplitudes. 

Studies on the matching between long- and short- distances in large N QCD, with 
the calculation of the Qj penguin matrix element and of at the NLO in the 1/N 
expansion have been presented in ref. [[37]]. However, a complete calculation of the 
K — > txtx matrix elements relevant to e'/e is not available yet. 

It has been recently emphasized |55| that cut-off based approaches should pay 



attention to higher- dimension operators which become relevant for matching scales 
below 2 GeV and may represent one of the largest sources of uncertainty in present 
calculations. On the other hand, the calculations based on dimensional regularization 
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Figure 3: Predicting e' je: a (Penguin) Comparative Anatomy of the Miinchen (dark gray) 
and Trieste (light gray) results in Fig. || (in units of 1CP 3 ). 



(for instance the Trieste one) may be safe in this respect if phenomenological input 
is used in order to encode in the hadronic matrix elements the physics at all scales. 

Lattice, as a regularization of QCD, is the first-principle approach to the prob- 
lem. Presently, very promising developments are being undertaken to circumvemt the 
technical and conceptual shortcomings related to the calculation of weak matrix ele- 
ments (for a recent survey see ref. |39fl). Among these are the Domain Wall Fermion 
approach [fH],41] which allows us to decouple the chiral symmetry from the continuum 
limit, and the possibility to circumvent the Maiani- Testa theorem [42| using the fact 



that lattice calculations are performed in finite volume |15fl . All these developments 
need a tremendous effort in machine power and in devising faster algorithms. Pre- 
liminary results for the calculations of both e'/e and the AI = 1/2 selection rule are 
expected during the next year. 



5 The AI =1/2 selection rule 

Without entering into the details of the various calculations I wish to illustrate 
with a simple exercise the crucial role of higher order chiral corrections (in general of 
non-factorizable contributions) for e'/e and the AI = 1/2 selection rule. In order to 
do that I focus on two semi-phenomenological approaches. 

A commonly used way of comparing the estimates of hadronic matrix elements 
in different approaches is via the so-called B factors which represent the ratio of 
the model matrix elements to the corresponding VSA values. However, care must be 
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Figure 4: Anatomy of the AI = 1/2 rule in the xQM. See the text for explanations. The 
cross-hairs indicate the experimental point. 



taken in the comparison of different models due to the scale dependence of the i?'s and 
the values used by different groups for the parameters that enter the VSA expressions. 
An alternative pictorial and synthetic way of analyzing different outcomes for e'/e is 
shown in Fig. |3], where a "comparative anatomy" of the early Trieste and Munchen 
predictions is presented. 

From the inspection of the various contributions it is apparent that the different 
outcome on the central value of e'/e is almost entirely due to the difference in the 
size of the Qq contribution. 

In the Munich approach [20] the AI =1/2 rule is used in order to determine 
phenomenologically the matrix elements of Qi,2 and, via operatorial relations, some 
of the matrix elements of the left-handed penguins. The approach does not allow for 
a phenomenological determination of the matrix elements of the penguin operators 
which are most relevant for e' / e , namely the gluonic penguin Q 6 and the electroweak 
penguin Qg. These matrix elements are taken around their leading 1/N values (fac- 
torization). 

In the semi-phenomenological approach of the Trieste group the size of the effects 
on the I = 0, 2 amplitudes is controlled by the phenomenological embedding of the 
AI =1/2 selection rule which determines the ranges of the model paremeters: the 
constituent quark mass, the quark and the gluon condensates. In terms of these 
parameters all matrix elements are computed. 

Fig. |] shows an anatomy of the xQM contributions which lead to the experimen- 
tal value of the AI = 1/2 selection rule for central values of the input parameters. 
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Point (1) represents the result obtained by neglecting QCD and taking the factorized 
matrix element for the tree-level operator Q 2 , which is the leading electroweak con- 
tribution. The ratio A /A 2 is thus found equal to y/2: by far off the experimental 
point (8). Step (2) includes the effects of perturbative QCD renormalization on the 
operators Qi,2 |H||- Step (3) shows the effect of including the gluonic penguin opera- 



tors ||45|| . Electroweak penguins |46| are numerically negligeable in the CP conserving 
amplitudes and are responsible for the very small shift in the A 2 direction. Therefore, 
perturbative QCD and factorization lead us from (1) to (4): a factor five away from 
the experimental ratio. 

Non-factorizable gluon-condensate corrections, a crucial model dependent effect 
entering at the leading order in the chiral expansion, produce a substantial reduction 



of the A 2 amplitude (5), as it was first observed by Pich and de Rafael p7| . Moving the 



analysis to 0(p 4 ), the chiral loop corrections, computed on the LO chiral lagrangian 
via dimensional regularization and minimal subtraction, lead us from (5) to (6), while 
the finite parts of the NLO counterterms calculated in the xQM approach lead to 
the point (7). Finally, step (8) represents the inclusion of ir-i]-i]' isospin breaking 
effects [0. 

This model dependent anatomy shows the relevance of non-factorizable contribu- 
tions and higher-order chiral corrections. The suggestion that chiral dynamics may 
be relevant to the understanding of the AI = 1/2 selection rule goes back to the 
work of Bardeen, Buras and Gerard |21| in the 1/N framework with a cutoff regu- 
larization. A pattern similar to that shown in Fig. ^ for the chiral loop corrections 
to A and A 2 was previously obtained, using dimensional regularization in a NLO 
chiral lagrangian analysis, by Missimer, Kambor and Wyler The Trieste group 
has extended their calculation to include the NLO contributions to the electroweak 
penguin matrix elements [ CT . 

Fig. |5] shows the contributions to e' / e of the various penguin operators, providing 
us with a finer anatomy of the NLO chiral corrections. It is clear that chiral-loop 
dynamics plays a subleading role in the electroweak penguin sector (Qs-io) while en- 
hancing by 60% the gluonic penguin (J = 0) matrix elements. The NLO enhancement 
of the Qq matrix element is what drives e'/e in the xQM to the 10~ 3 ballpark. 

As a consequence, the xQM analysis shows that the same dynamics that is relevant 
to the reproduction of the CP conserving A amplitude (Fig. |J) is at work in the CP 
violating sector, albeit with a reduced strenght. 

In order to ascertain whether these model features represent real QCD effects 



we must wait for future improvements in lattice calculations |p9|| . On the other 
hand, indications for such a dynamics follow also from the analitic properties of the 
K —>■ 7T7T amplitudes ||. It is important to notice however that the size of the effect 
so derived is generally not enough to fully account for the AI = 1/2 rule. Other non- 
factorizable contributions are needed to further enhance the CP conserving 1 = 
amplitude, and to reduce the large 1 = 2 amplitude obtained from perturbative QCD 
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Figure 5: Anatomy of e'/e (in units of 10 -3 ) in the xQM approach. In black the LO 
results (which includes the non-factorizable gluonic corrections), in half-tone the effect of 
the inclusion of chiral-loop corrections and in light gray the complete 0(p 4 ) estimate. 



and factorization. In the xQM approach, for instance, the fit of the AI = 1/2 rule 
is due to the interplay of NLO chiral corrections and non-factorizable soft gluonic 
contributions (at LO in the chiral expansion). 

6 Conclusions 

In summary, those semi-phenomenological approaches which reproduce the AI = 
1/2 selection rule in K — > tttt decays, generally agree in the pattern and size of the 
1 = 2 hadronic matrix elements with the existing lattice calculations. On the other 
hand, the A J =1/2 rule forces upon us large deviations from the naive factorization 
for the 1 = amplitudes: B— factors of O(10) are required for (<5i,2) - Here is were 
lattice calculations presently suffer from large sistematic uncertainties. 

In the Trieste and Dortmund calculations, which reproduce the CP conserving 
K —>■ tttt amplitudes, non-factorizable effects (mainly due to final-state interactions) 
enhance the hadronic matrix elements of the gluonic penguins, and give B e /B^ « 2. 
Similar indications stem from recent 1/N |33| and dispersive || approaches. The 
direct calculation of K — > tttt amplitudes and unquenching are needed in the lattice 
calculations in order to account for final state interactions. Promising and exiciting 
work in this direction is in progress. 
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